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In the VKS2 (von Karman Sodium 2) successful dynamo experiment of September 2006, the 
magnetic field that was observed showed a strong axisymmetric component, implying that non 
axisymmetric components of the flow field were acting. By modeling the induction effect of the 
spiraling flow between the blades of the impellers in a kinematic dynamo code, we find that the 
axisymmetric magnetic mode is excited and becomes dominant in the vicinity of the dynamo thresh- 
old. The control parameters are the magnetic Reynolds number of the mean flow, the coefficient 
measuring the induction effect, a, and the type of boundary conditions (vacuum for steel impellers 
and normal field for soft iron impellers). We show that using realistic values of a, the observed 
critical magnetic Reynolds number, Rm^ ^ 32, can be reached easily with ferromagnetic boundary 
conditions. We conjecture that the dynamo action achieved in this experiment may not be related 
to the turbulence in the bulk of the fiow, but rather to the alpha effect induced by the impellers. 

PACS numbers: 47.65.-d, 52.65Kj, 91.25Cw 



The interest of the scientific community for the dy- 
namo action in liquid metal has been renewed since 2000 
in the wake of successful experiments 0, We fo- 
cus in this Letter on the Cadarache VKS2 (von Karman 
Sodium 2) experiment Q where the dynamo effect oc- 
curred beyond the critical magnetic Reynolds number 
Rm^ ^ 32. The growing magnetic field that was observed 
above threshold was mainly axisymmetric. Dynamo ac- 
tion was found with soft iron impellers but did not occur 
with stainless steel impellers, using the same available 
power. The purpose of the present Letter is to present 
arguments based on modeling and numerical simulations 
to justify the observed axisymmetric mode. The two key 
ingredients are the so-called a-effect and ferromagnetic 
boundary conditions. 

From Cowling's theorem 0] the axisymmetric part of 
the fiow cannot be responsible alone for the generation of 
an axisymmetric field. In the experiment, the deviation 
from axisymmetry of the velocity field is thus expected 
to participate to the dynamo process, generating an ax- 
isymmetric electromotive force somewhere in the fiuid 
and triggering the dynamo growth. Two other magneto- 
hydrodynamic experiments using counter-rotating pro- 
pellers have also led to this conclusion either on observa- 
tional^, Q or on numerical grounds |7|. It was argued 
in [l, 0^0 that turbulence effects were responsible for 
the occurrence of the axisymmetric magnetic field. In Q 
the turbulence effects were held responsible for the dy- 
namo threshold in VKS2. In addition an a-effect close to 
the discs was invoked to justify the axisymmetry of the 
generated field. 



The VKS2 configuration is modeled as follows: The 
computational domain is divided into three concentric 
cylinders of axis embedded in an insulating sphere 
(figure [T]). The inner cylinder r < Rq contains the mov- 
ing sodium of conductivity ctq. The cylindrical shell 
Rq ^ r < Ri contains stagnant sodium (ai = ctq). 
The cylindrical shell Ri < r < R2 is a. layer of cop- 
per ((72 = 4.5<Jo). We choose Rq as unit length and the 
geometric parameters are Ri/Rq = 1.4, R2/R0 = 1-6 
and Rout/Ro = IO7 with the same height H/Rq = 1.8 
for the three cylinders. The velocity field of the liquid 
sodium U {r < Rq) is axisymmetric and is obtained from 
time averaged measurements of a water experiment in 
the VKS2 setting The meridional fiow field and the 
isolines of the azimuthal component of the velocity are 
shown in figure [H The fiow field is maintained by two 
counter-rotating impellers located at both ends of the 
vessel. They each occupy a cylindrical volume of radius 
Ro and thickness H/20. 

The only non axisymmetric parts of the device are the 
blades composing the impellers (eight each), and we rep- 
resent the effect of these blades by adding the electromo- 
tive force a{'B-e0)e0 to the induction equation. This force 
is assumed to act only in the two fiuid cylinders occupied 
by the blades as they rotate. This model is deduced as 
follows: As the blades rotate, the fiuid is ejected radi- 
ally outwards and radial vortices are created due to the 
rotation rate gradient, see the white arrows in figure O 
We then assume that when a magnetic field is applied, 
the mean induced current is dominated by its azimuthal 
component and that the shearing action of the vortices 



FIG. 1: Geometry of the computational domain in a merid- 
ional plane. The axisymmetric meridional flow {Ur,Uz) is 
represented on the left with arrows. On the right the isoval- 
ues oi Ue are plotted. The radii Ri and height H as well as 
the conductivities ai are defined in the text. 

is cumulative. The axial current is expected to vanish in 
the vicinity of the discs and, as a result, we set the azz 
component of the a tensor to be zero. This phenomenon 
is illustrated on fig. [2l As the magnetic lines (B • eo)e0 
(thin line in figure [2]) are distorted by the helical vor- 
tices, small scale magnetic loops are created, which in 
turn create an azimuthal electromotive force a^B ■ eo)eo^ 
where the parameter a is negative and depends on the 
impeller rotation rate (see below for an estimate), the 
size, curvature and number of blades. We thus attempt 
to represent the complexity of the induction effects in a 
poorly known flow close to a boundary by a single scalar 
parameter. The handedness of the helical vortices gen- 
erated between the blades is directly linked to the rota- 
tion sign of the impellers: orienting the rotation axis of 
the bottom impeller towards the center of the container, 
a positive rotation drives a set of eight radial positive 
helices creating a negative a coefficient. The counter- 
rotation configuration with (^bot = — ^top is invariant by 
rotation of tt about any horizontal axis [10], therefore 
the two impellers produce the same a. Actually, what- 
ever the sign for the rotation of the bottom impeller, 
the product Qa is always negative. Henceforth we set 
(^bot = > 0. 

We follow the standard kinematic dynamo approach 
by solving the induction equation 

dtB = V X (UxB + a(B-e^)e^)- Vx(7^VxB). (1) 

where 77 = 1/cr/io, Mo being the vacuum magnetic per- 
meability, a the fluid conductivity and a the unknown 
parameter which models the induction effect in the vo^j 
ume occupied in average by the impellers. The field U 
is equal to U in the cylinder r < Rq and is zero in the 
two cylindrical shells Rq < r < Ri^ Ri < r < R2. The 
conductivity field is such that a = ctq in the inner cylin- 
der r < Ro and in the stagnant sodium shell, and a = a2 
in the copper layer. The conductivity jump at r = i?i 
is accounted for by enforcing the continuity of the tan- 
gent component of the electric field. Note that the flow 



FIG. 2: Sketch of the top impeller with eight radial blades 
attached to the disc. Due to the rotation of the disc (gray 
arrow on the top) , radial helical vortices are produced between 
the blades (thick white arrows). An axisymmetric azimuthal 
field {'B-ee)ee (thin Une) is distorted leading to eight magnetic 
loops. 

behind the impellers has not been simulated (see dis- 
cussion below). To simplify the parameterization of the 
alpha effect, a vanishes outside two thin cylindrical lay- 
ers of radius Rq and thickness H/ 20 situated at each end 
of the container and its intensity varies smoothly within 
these layers using a tanh regularization function. The 
magnetic Reynolds number is Rm = cfqiiqUqRq where 
[/o is the maximum speed of the flow U. 

Two types of boundary conditions are compared 
to model the experiment: they are referred to as 
insulating (I) and ferromagnetic (F) boundary con- 
ditions. We say that insulating boundary con- 
ditions (I) are used when the continuity of the 
magnetic field is enforced at the interface between 
the conducting regions and the vacuum (i.e., at 
\z\ = H/2, 0<r <R2 and at < \z\ < H/2, r = R2). 
We say that ferromagnetic boundary conditions (F) 
are used when we enforce B x = at 

\z\= i^/2, < r < and we require B to be contin- 
uous across the other boundaries < \z\ < i^/2, r = R2 
and 1^1 = H/2, Rq <r < R2. The condition B x = is 
meant to mimic the effect of ferromagnetic discs with infi- 
nite magnetic permeability located at the top and bottom 
of the vessel [jjj . These boundary conditions model steel 
impellers (I) and soft iron impellers (F), respectively. 

Our numerical code [I2I, [13] uses a finite element 
Galerkin method for the spatial discretization in the 
meridional plane (r, z) and a Fourier series decomposi- 
tion along the azimuthal direction. Since the coefficients 
of the induction equation do not depend on the azimuthal 
angle, the different azimuthal modes (m = 0,1,2,...) 
are decoupled. The time is discretized with a semi- 
implicit backward finite difference method of second or- 
der. Growth rates are computed for the two types of 
magnetic boundary conditions, (I) or (F), and for vari- 
ous values of a and Rm. The zero growth rate condition 
defines a critical curve in the (a, Rm) plane. A compila- 
tion of results obtained for the (I) and (F) cases and for 
Fourier modes m = and m = 1 is shown in figure [3l 
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FIG. 3: Dynamo threshold Rrrf versus —a, for (x, A) m = 0, 
m = 1 and where (F) and (I) indicate ferromagnetic and 
insulating boundary conditions respectively. 




FIG. 4: Top: Geometry of the m = (F) growing magnetic 
field for i?m^ = 32, |a| = 2.2 > |ac| = 2.1 . The three 
cylinders are also represented. The lines correspond to the 
positive isovalues of Be and the arrows to the meridional field 
(Br,Bz). Bottom: isovalue of 25% of the maximum energy 
plotted in the 3D domain, together with some field lines. 



In order to make comparisons with experimental re- 
sults published in we first focus our attention on the 
axisymmetric magnetic mode (m = 0), which was a pri- 
ori eliminated in former studies by Cowling's theorem. 
The two critical curves corresponding to the boundary 
conditions (I) and (F) and the Fourier mode m = 
reveal that the dynamo action occurs only for a neg- 
ative, thus confirming the phenomenological argument 
discussed above to justify the sign of a. Both curves 
clearly show that increasing \a\ lowers the critical mag- 
netic Reynolds number. Moreover, the (F) boundary 
condition yields significantly smaller thresholds. This ob- 
servation agrees with previous results obtained in a dif- 
ferent context [3, [m . Both curves can be represented 
by the scaling law Rm'^ ~ Ici^l"'^? where q ~ 0.67, sug- 
gesting that the dynamo process is similar for both types 
of boundary conditions. 

The value of a corresponding to the dynamo threshold 
in VKS2, i?m^ = 32, is found to be |a| = 2.1 for the (F) 
boundary condition and |a| = 3.3 for the (I) boundary 
condition. Since the steel impellers have not produced 
dynamo action, this suggests that, as a first guess, the 
effective |a| for VKS2 is between 2.1 and 3.3. A rough 
estimate of a is given by a ~ u^h/rj where is a typ- 
ical flow intensity between the blades and h = H/20. 
From we estimate u/Uq ^ 0.2 at the impeller half ra- 
dius. Then, for Rm = 32, 7^=0.1 m^.s~^ and Rq = 0.2m, 
we find a ~ 1.8m.s~^, making the value \a\ = 2.1m.s~^ 
plausible. If |a| is indeed close to 2.1, the critical value of 
Rm for the (I) boundary condition is close to 43, which 
apparently contradicts the experiments since Rm = 50 
has been reached without dynamo action using steel im- 
pellers. To explain this discrepancy we recall that it has 
been shown in fid\ that for the Fourier mode m = 1 the 
flow of liquid sodium behind the impellers acts against 
the dynamo action. We expect this anti-dynamo effect to 
be active also for the Fourier mode m = 0, but we did not 
include this extra layer of fluid in the present study. The 
negative consequences of this secondary flow are proba- 
bly screened in the case of soft iron impellers. The nu- 
merical evidence justifying this assertion has to await the 
availability of a code describing conducting domains with 
different magnetic permeabilities (work in progress). 

The stationary eigenvector corresponding to the (F) 
boundary condition for m = and \a\ = 2.2 is shown 
in figure [H The magnetic field is mainly concentrated 
along the z-axis in the fluid region r < Rq^ and it is 
azimuthally dominated in the plane z = for r > Rq. 
The radial component is odd with respect to z whereas 
the azimuthal and vertical components are even and of 
opposite sign. These features are compatible with the 
magnetic field measured at saturation in the experimen- 
tal dynamo regime obtained with soft iron impellers 
The eigenvector corresponding to the (I) boundary con- 
dition for m = and \a\ = 2.2 (not shown here) is as- 
sociated to an eigenvalue with a nonzero imaginary part. 
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The resulting dynamo is periodic with a reversing dipo- 
lar moment. This result, not detailed here, is in strong 
contrast with the case a = (I), for which a non os- 
cillatory transverse dipole mode (m = 1) with internal 
banana-like structures is the only unstable mode 0. 

We have also investigated the effect of the a model 
on the Fourier mode m = 1. This mode has not been 
observed in the VKS2 experiment although it has been 
predicted to be the most unstable in a xisy mmetric kine- 
matic dynamo simulations [1, [H, [13, [3 • figure [31 
for the mode m = 1, we see that the critical magnetic 
Reynolds number grows sharply with |a| for a < and 
decreases with \a\ for a > 0. We are thus led to con- 
clude that the growth of the m = 1 mode is hindered 
by the very induction effect that generates the m = 
mode! This surprising result may explain the absence of 
the m = 1 mode in the experiment, at least for Rm < 50. 

The present study shows that the positive results of the 
dynamo experiment of September 2006 (VKS2) as well 
as the former negative results (VKSl) may be explained 
by a simple mean induction effect induced by the radial 
vortices trapped between the blades of the impellers. We 
conjecture that turbulence may not play a role as essen- 
tial as initially believed in the VKS experiments. 

This work also shows the importance of using ferro- 
magnetic materials, corresponding to the (F) boundary 
conditions in fig.[3l A research program aiming at explor- 
ing the saturation regime using nonlinear computations 
and materials with heterogeneous magnetic permeabili- 
ties is currently engaged. 
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